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Abstract 

We review the theory of hadronic atoms in QCD+QED. The non-relativistic effective Lagrangian 
approach, used to describe this type of bound states, is illustrated with the case of tt^-k^ 
atoms. In addition, we discuss the evaluation of isospin-breaking corrections to hadronic atom 
observables by invoking chiral perturbation theory. 
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Hadronic atoms are bound states of hadrons, held together predommately by 
the static Coulomb force. Simple examples are pionium, a bound state of two 
pions with opposite electric charge (7r+,7r~), and pionic hydrogen, a bound state 
of a 7r~ and a proton. Pionium is the analogue of positronium in Quantum 
Electrodynamics. A more complex example is pionic deuterium ~ a Coulombic 
bound state of a 7r~ and a deuteron. The latter itself is a composite state of 
a proton and a neutron, bound at much smaller distances than the size of the 
hadronic atom. 

The average distance between the constituents of a hadronic atom is given by 
the Bohr radius 

1 M1M2 , , 

a/Xc Ml + M2 

where a ~ 1/137 is the fine-structure constant, fj,c stands for the reduced mass, 
and Ml, M2 are the masses of the constituents. Typically, the Bohr radius is of the 
order of a few hundred Fermi, much larger than the range of strong interactions. 
Individual hadrons in the atom spend most of the time at distances where strong 
interactions are practically absent. For this reason, observables of hadronic atoms 
are barely affected by the strong interactions. 

If the interaction between the constituents were purely electromagnetic and 
non-relativistic, the energy levels of the atom would be given by the standard 
quantum-mechanical formula 

En = Ml + M2 - ^ fica\ n = l,2,---. (1.2) 

Aside from relativistic corrections which generate higher order terms in a, this 
formula is modified in the presence of strong interactions in two ways. First, the 
energy levels are shifted from their purely electromagnetic value. Furthermore, 
because the atoms can decay also via strong interactions (example: the decay of 
pionium into a neutral pion pair through the charge-exchange reaction 7r~^7r~ — > 
), the energy levels are broadened. The effect on the ground-state energy 
level is illustrated in Fig. [H The pionium lifetime in the ground state, r = l/F ~ 
3 X 10^^^ s, is still much smaller than the charged pion lifetime — 10~^ s. 
Despite the short lifetime of the atom, the pions travel many times around each 
other before the atom decays, as the ratio ^ficO^/T ~ 8 x 10'^ indicates. As 
a consequence of this, pionium can be considered a quasi-stable bound state 
with a clearly defined structure of (almost Coulombic) energy levels. The same 
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Figure 1: Schematic representation of the shift of the ground state energy level in 
pionium. AE^^^ and T denote the strong energy shift and the width in the ground 
state (we omit indices in F throughout, because we will consider the widths of 
ground states only) . The symbols em and meas denote the purely electromagnetic 
and the measured energy levels, respectively. 

statement is valid for many other hadronic atoms. 

Because the size of hadronic atoms is much larger than the range of strong 
interactions, the energy levels of the atoms can depend only on the characteristics 
of hadronic interactions at asymptotically large distances. These are usually 
described in terms of the parameters in the effective range expansion: scattering 
length, effective range, shape parameters. The situation is analogous to the 
calculation of the classical static electric field generated by a charge distribution: 
at asymptotic distances, the electric field depends only on the multipole moments 
which describe the charge distribution. 

Deser, Goldberger, Baumann and Thirring (DGBT) were the first to derive - 
at leading order in the fine-structure constant a - a formula for the complex shift 
of the energy level of a hadronic atom The real and imaginary parts of this 
shift define the displacement and the width of a given level, generated by the 
strong interactions. The formula for the ground state reads 



where T denotes the complex elastic scattering amplitude of the constituents at 
the threshold. The ellipses stand for higher order isospin breaking corrections, 
which will be discussed in detail later in this article. The formula can be trivially 
generalized to the case of excited energy levels. 

The DGBT formula (jl.3p plays a central role in the theory of hadronic atoms, 
because it allows one to extract the threshold amplitude T from the experimen- 
tally measured energy and width of the atom. Further, the real and imaginary 



AE""^' - - F 



2a^fil T + ■■■ 



(1.3) 
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parts of T are related to the hadronic scattering lengths. From this we conclude 
that the experimental study of hadronic atoms provides us with a source for the 
determination of these scattering lengths. The above-mentioned huge difference 
in the atomic and strong interaction scales is very advantageous in the present 
context, since the atomic observables depend (at leading order) exactly on those 
quantities (scattering lengths) which one wants to extract from the experiment - 
they are not sensitive to the short-range details of strong interactions. 

In most cases, the first term on the right-hand side of Eq. (jl.3p does not match 
the available experimental precision - next-to-leading order corrections, indicated 
with the ellipses in Eq. (jl.3p , are needed as well for a precise determination of the 
scattering lengths. The aim of any theory of hadronic atoms must be to provide 
a systematic framework for the calculation of these corrections. Here, we will 
carry out the calculations by using a non-relativistic effective theory of Quantum 
Chromodynamics -|- Quantum Electrodynamics (QCD+QED). We will illustrate 
the method by means of pionium decay, and will briefly consider the application 
of the same approach to pionic hydrogen and pionic deuterium. We no not discuss 
quantum-mechanical potential models, because these methods introduce inherent 
model-dependent artefacts which cannot be controlled. 

The article is organized as follows. In section [2] we consider the physics back- 
ground behind the experiments on various hadronic atoms. In section[3]we briefly 
discuss the essentials of the non-relativistic effective Lagrangian approach, which 
is our tool to describe hadronic atoms. Section H] forms the backbone of the 
present article. In this section we construct, step by step, the effective field the- 
ory approach to pionium decays. The same approach is applied in section [5] to 
the description of pionic hydrogen and of pionic deuterium. Finally, section [6] 
contains a brief summary and outlook for future research in the field. 



2 Physics background 



Several experiments with hadronic atoms are presently running. The DIRAC 
collaboration at CERN is measuring the lifetime of pionium (3, 3, Q, Q, Q, Q, S) 
and plans to determine the lifetime of irK atoms as well (^). The Pionic Hydrogen 
collaboration at PSI (|lol. [ill, [l^. [l^. 1^ studies the spectrum of pionic hydrogen 
and pionic deuterium, whereas the DEAR/SIDDHARTA collaboration at LNF- 
INFN (Q, Q, 0, Q) plans to determine the ground state energy and width of 
kaonic hydrogen at a much better accuracy than in previous experiments carried 
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out at KEK (Q,Q). In addition, SIDDHARTA plans the first ever measurement 
of the spectrum of kaonic deuterium. 

These experiments eventuahy result in a precise determination of various hadro- 
nic scattering lengths. Let us recall why the results will be important for the 
investigation of several fundamental properties of QCD. 

1. We start with the DIRAC experiment at CERN fl Q, Q, 3, Q, 3)- The 
decay width of the ground state of pionium into a vr'^vr*' pair is related to the 
difference of the S-wave vrvr scattering lengths ao,a2 with total isospin and 2, 

„ 2 



aV(ao - 02)^ + • 



9 



(2.1) 



Here, p* = (M^ — M^q — j M^a^)^/^H is the CM momentum of the neutral pion 

pair after decay, Mt^jM^o are the charged and neutral pion masses, respectively, 
and the ellipses stand for terms of higher order in isospin breaking. 

It is expected that the DIRAC experiment will finally provide a value for 
|ao — which is accurate up to a few percent. Other experiments, where the vrvr 

or 
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scattering lengths are determined from K^i decays (21. H, [23I, 24 . [2^^ or from 
studying the cusp structures in K ^ 3tt decays (0, 
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31 



32, 



33, 



35I ). yield competitive results in terms of accuracy. On the other hand, the 



difference an — go is particularly sensitive to the value of the quark condensate 



in QCD 136 
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33 ). In the so-called "standard" scenario which assumes a 



large condensate, the expansion of the pion mass in terms of the quark mass is 



Mi = 



327r2F2 



m 



{rriu + md), (2.2) 



where the second term on the right hand side of the first equation is small MQ ). 
Here, F is the pion decay constant F^^ in the chiral limit, mu,md are the light 
quark masses. Is denotes one of the low-energy constants (LECs) in chiral per- 

and the quantity B is 



turbation theory (ChPT, see, e.g.. References ([41 



related to the quark condensate in the chiral limit ([40|, |42| ). Further, if chiral 



symmetry breaking in QCD proceeds according to the standard picture, a very 



accurate description of the scattering lengt 
2-loop ChPT with the Roy equations (l43l . 



IS ao , a2 can be achieved by combining 



4J) 



ao = 0.220 ± 0.005 , 02 = -0.0444 ± 0.0010 , ao - a2 = 0.265 ± 0.004 . (2.3) 



Equipped with this precise theoretical prediction, one may perform a direct ex- 
perimental test of the chiral symmetry breaking scenario in QCD. Namely, if the 
measured value of ao — a2 differs significantly from the theoretical prediction. 
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this would suggest that chiral symmetry breaking proceeds in a manner which is 
different from the standard scenario. The present situation concerning the verifi- 
cation of the predictions (j2.3p is the fohowing. Lattice results for 02 (|45|, 1461 ) agree 
with the prediction within one standard deviation, see Ref. (j46l ) for a compilation 
of predictions, lattice calculations and data. Due to technical difficulties (discon- 
nected graphs), oq has not yet been measured with this technique. On the other 
hand, using the LECs ^3^4 (or their SU (3) x SU{3) counterparts) determined from 
the lattice and converting these into a value of gq again leads to agreement with 



the prediction (j2.3p , within one standard deviation (j47|, |4^, |49|, |50|, |51|, |5^, l53l ) . 
We refer the interested reader to Ref. (|48l ) for discussions and for a review. On 
the experimental side, data of the DIRAC collaboration on pionium lifetime 



of NA48/2 on the cusp in — > Svr decays (l29l ) and on K^a events (l23l ) (applying 
isospin-breaking corrections as described in Refs. (12J, |25|)) neatly confirm the 
predictions, although partly still with considerable uncertainties. 

2. Next, we briefly consider the proposed measurement of the ttK atom life- 

one to extract the value of the isospin-odd S-wave ttK 
scattering length Oq . The calculations of this scattering length, carried out in 
ChPT up to two loops, lead to a rather contradictory picture: it turns out (|54l .l55l) 
that the two-loop contribution to this quantity is apparently larger than the one- 
loop correction. On the other hand, the result at two loops agrees with the 
analysis carried out on the basis of Roy equations (jSGl ) . The situation is puzzling, 
because, if correct, the convergence of the ChPT series for pion-kaon scattering is 
under question. It is clear that a precise knowledge of the experimental value of 
the scattering length is an important ingredient to the solution of this puzzle. For 
more comments concerning this point, we refer the interested reader to section 2 
of the review (1571 ). 

3. From the measurement of the pionic hydrogen erier gy s hift and width by 
the Pionic Hydrogen collaboration at PSI (jld . 
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12 



13l . Il4l ). one can extract 
the isospin even and odd S-wave vrA^ scattering lengths a^^ and a^_|_. Using 
Effective Field Theory methods (EFT) in the two-nucleon sector, one can also 



relate the pion-deuteron scattering l eng ths to t 



61 



63 



64], 



65 



67 



68|, 



re pion-nucleon ones. [See, e.g.. 



A similar result 



69 



0, 0, 0; 



Refs. (158|, 159|, 1& 

can be obtained with quantum-mechanical multiple-scattering theory (see, e.g., 
Refs. (0, 0, 0)).] Thus, the measurement of the energy shift and width of 
pionic deuterium results in additional constraints on the values of aQ_^ and Oq^. 
ttN scattering lengths are quantities of fundamental importance in low-energy 
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hadronic physics by themselves, since they test the exact pattern of exphcit chiral 
symmetry breaking. Moreover, knowledge of the exact values of the scattering 
lengths also affects our understanding of more complicated systems where the 
vrA^ interaction serves as an input, e.g. NN interaction, pion-nucleus scattering, 
three-nucleon forces, etc. In addition, high-precision values of the ttN scattering 
lengths are used as an input for the determination of different basic parameters 
of QCD at low energies more accurately. One example is the ttNN coupling 
constant gTrN iy, which is obtained from the Goldberger-Myazawa-Oehme (GMO) 



sum rule (76 



77 



78l ) , where a particular combination of scattering lengths enters 



as a subtraction constant. Other important quantities, which can be obtained 
by using the 5- wave ttN scattering lengths as an input, are the so-called pion- 
nucleon sigma-term and the strangeness content of the nucleon. The sigma-term 
o'ttN, which measures the explicit breaking of chiral symmetry in the one-nucleon 
sector, is defined by 

cr-n-Af = —^—(ps\m(uu + dd)\ps) , (2.4) 

where \ps) denotes a one-nucleon state, with momentum p and spin s, and tun is 
the nucleon mass. The sigma-term is related to the strangeness content y of the 
nucleon, and to the SU{3) symmetry breaking part of the strong Hamiltonian, 

"^■^ {ps\uu + dd — 2ss\ps) = ( — 1 ) (1 — y)cFTrN , 



2mN V 

y - ^^^^^ 

{ps\uu + dd\ps) ' 

where m<j denotes the strange quark mass. In the analysis of the experimental 
data, one uses 5-wave vrA^ scattering lengths as input in the dispersion relations, 
which provide the extrapolation of the isospin even pion-nucleon scattering am- 
plitude from threshold down to the Cheng-Dashen point. We refer the interested 

n 

reader to Ref. (j79l ) for details. In this reference, the value u-^-n — 45 MeV was 
obtained. [In Ref. (jSOl ). a value for the sigma-term which is considerably larger 
than 45 MeV is claimed to follow from more recent data.] The sigma-term is 
rather sensitive to the scattering lengths Consequently, an accurate mea- 

surement of the latter will have a large impact on the experimentally determined 
values of a-itN and y. Finally, we note that the sigma-term is accessible through 
lattice calculations, see e.g. Ref. (ISlI ) and references cited there. It even plays a 
role in astrophysical applications. As an example for such an impact, we refer the 
interested reader to the recent publication (|82l ) and the references given there. 
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4. Last but not least, we discuss the DEAR/SIDDHARTA experiment at LNF- 
INFN (|l3,Q,0,Q). It plans to determine KN scattering lengths from data 
on kaonic hydrogen and on kaonic deuterium atoms. We believe that it would 
be very useful to carry out a comparison of the scattering lengths so determined 
with different theoretica.1 predictions based on the unitarization of the lowest 
order ChPT amplitude dssl, 
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90 



91 



Indeed, it turns out 



that even the data from kaonic hydrogen alone impose rather stringent constraints 
on the values of the KN scattering lengths. In some cases, DEAR/SIDDHARTA 
data seem not to be compatible with the scattering sector (StI. [9I). It is clear 
that imposing additional constraints from Kd data makes the issue even more 
pronounced. In our opinion, it is important to check whether the unitarization 
approach passes this test. 



3 The non-relativistic effective theory 

At leading order, the DGBT formula in Eq. (jl.3p is universal: it looks exactly the 
same in potential scattering theory (where it was derived first) and in quantum 
field theory. This fact is due to the huge difference between the atomic and the 
strong interaction scales mentioned in the introduction. On the other hand, the 
isospin-breaking corrections to this relation, which are due to electromagnetic 
interactions and to the quark mass difference — are not universal. In this 
article, we describe a systematic theory of hadronic atoms within QCD+QED, 
which enables one to calculate these corrections in a simple and elegant manner, 
with an accuracy that matches the experimental precision. [ Because ChPT is the 
low-energy effective theory of QCD+QED, one might be tempted to start from 
this framework. However, describing bound states in ChPT by using standard 
techniques, based on the Bethe-Salpeter equation or on 3-dimensional reductions 



thereof ((92, 



93 



94 l95l . l96l . l97l ). is a complicated enterprise, which makes it very 
difficult to reach the required precision. We do not, therefore, discuss this method 
here.] 

Our framework is based on the existence of several different momentum scales in 
the problem. Counting powers of the fine-structure constant a, we have to assign 
the order a" to the scale related to the pion mass, because has a non vanishing 
value also in the absence of electromagnetic interactions. On the other hand, the 
momentum scale corresponding to atomic phenomena is given by the inverse Bohr 
radius - i.e., the average 3-momenta inside the atom are pav — ^ij^ = a/^'o and 
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count as order a. From this one concludes that a non-relativistic approach, based 
on an expansion in (small) momenta, is the appropriate framework to describe 
hadronic atoms, because the momentum expansion translates into an expansion 
in the fine-structure constant for hadronic atom observables. The advantage of 
considering a non-relativistic framework consists in the simple treatment of bound 
states: they can be described by the Schrodinger equatiori^. 
Let us list some very general properties of this approach. 

i) The framework uses the language and methods of (non-relativistic) quan- 
tum field theory. In particular, the calculations are based on effective La- 
grangians and Hamiltonians. 

ii) The non-relativistic approach allows one to keep the number of heavy par- 
ticles conserved, by construction. In other words, one always stays within 
a restricted sector in Fock space. 

iii) The non-relativistic theory describes matrix elements at small external mo- 
menta. All high-energy effects ~ like transitions to sectors with a different 
number of heavy particles - are encoded in the coupling constants of the 
effective Lagrangian, which are determined through matching to the un- 
derlying theory. In this manner, one makes sure that the effective and the 
underlying theory are equivalent at low energies. 

iv) Power counting rules are at the heart of any effective field theory. The non- 
relativistic power counting at tree-level amounts to counting the number of 
space derivatives in various terms. Because each non-relativistic momentum 
is of order of a/Xc, the contributions to the bound-state energy from terms 
containing higher derivatives are suppressed by additional powers of a. To 
carry out calculations of the bound-state energy spectrum at a fixed order in 
a, a finite number of terms in the Lagrangian thus suffices [for comparison, 
in ChPT the number of the relevant terms is infinite]. 

v) The non-relativistic Lagrangian is used to generate Feynman graphs in a 
standard manner. Strong loops respect the power counting, if dimensional 
regularization is used. Loops with photons can also be made consis tent with 
power counting by applying the so-called threshold expansion (9^, 100 ). 

vi) It is useful to extend the power counting to include the isospin break- 



^ Caswell and Lepage (l98f ) were the first to use a systematics non-relativistic effective La- 
grangian approach to investigate bound states in QED. 
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ing effects which are generated by the quark mass difference — m„ as 
well - along with the electromagnetic corrections characterized by the fine- 
structure constant a. There is no strict rule for doing this. We now note 
that the effect of md — mu in the pion mass is of order {md — mu)'^, and linear 
for kaons and nucleons. It is therefore convenient to introduce a common 
isospin-breaking parameter 5 and count a ~ (m^ — m^j)^ ~ 5 in pionium, 
~ (^d ~ iT^u) ~ S otherwise. This has the advantage that the leading 
corrections to the hadronic atom observables, generated by a and md — rriu, 
are counted at the same order in 5. 

vii)) At the end of the day, when the hadronic atom spectrum is calculated and 
the matching to the underlying theory is performed, there is no reference 
left to the non-relativistic theory in the final result. The non-relativistic 
approach is used only at an intermediate stage, in order to facilitate the 
calculations. 

Our main goal here is to first evaluate the 0{5) isospin-breaking corrections to 
the leading order strong energy shift and width. These corrections are indicated 
by the ellipses in the DGBT formula Eq. (jl.3p . Due to lack of space, we con- 
centrate on those corrections that are relevant for the width of the ground state. 
These are more easy to pin down than those for the real part of the energy shift. 
In a second step, the right hand side of equation Equation (|1.3p will be expressed 
in terms of isospin symmetric scattering lengths, up to isospin breaking correc- 
tions. The path to a comparison of the DGBT formula with experimental data is 
then paved, and a precise determination of scattering lengths becomes feasible. 

In the next section, we consider in some detail the construction of a non- 
relativistic theory along these lines. The framew o rk was d e veloped during the 




last decade in Ref s . (157 
110 . 111 . 112 . II3I . 114 . 
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4 Pionium: Decay of the ground state 



Instead of presenting the non-relativistic effective theory in its full generality, 
we have decided to explain the method with one particular example, the decay 
of pioniu m. Tech nical details will be skipped - these can be found, e.g., in 
Refs. (|57l . ll04l . ll09l ). For a thorough discu ssion of the properties of non-relativistic 
theories, we refer the reader to Refs. (|57l . Ilia ). 
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Figure 2: The singularity structure of the tttt partial- wave scattering amplitudes 
in the complex s-plane. The shaded area denotes the low-energy domain. 

4.1 Non-relativistic framework: strong sector 

We start with a non-relativistic theory for pions, in the absence of photons, which 
will be included afterwards. On the other hand, it is very convenient to keep from 
the beginning the masses of charged and neutral pions at their physical values. 
This is a perfectly consistent procedure, because in the non-relativistic theory, 
these masses are not renormalized, even when the electromagnetic interactions 
are turned on. 

Our starting point is the relativistic amplitude for the process vr"(pi)7r''(p2) 
tt'^ {ps)tt'^ (pi) , where a,b,c,d = it or 0. Performing a partial-wave expansion, 
we arrive at the partial-wave amplitudes that depend on a single Mandelstam 
variable s = {pi +P2)^- Assuming further s to be a complex variable, we consider 
the singularity structure of the partial-wave amplitudes in the low-energy region 
\s — 4M^| <C [For definiteness, we consider here the sector with total charge 
Q = 0. Other sectors can be discussed analogously.] As it is well known, the 
partial-wave amplitude is holomorphic in the complex s-plane, cut along the 
positive real axis for s > 4M^o > see Fig. [2j Another branch point corresponding 
to the two charged pion threshold is located at s = 4M^, whereas the first inelastic 
threshold is located at s = 16M^o • In addition, there is a cut on the negative real 
axis. However, the distance between these faraway singularities and the 2-pion 
threshold is of the order of the pion mass squared. Consequently, in the low- 
energy region, which includes the neutral and charged two-pion thresholds (far 
below the first inelastic threshold), the partial- wave amplitude has a particularly 



simple form (j24l ). 



Ti{s) = Ms) + iBi{s)a{s) + iCi{s)aois) + Di{s)a{s)ao{s) , 



a{s) = yi-^, ao(.) = yi-^, (4.1) 
where Ai{s), - ■ ■ Di{s) are meromorphic functions in the low-energy domain. 
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We now present a framework which describes the relativistic vrvr scattering am- 
pUtude in the low-energy region, and thus reproduces this structure of the am- 
pUtude. The kinetic term in the Lagrangian is fixed through the non-relativistic 
expansion of the relativistic one-particle energy, 

where ^*±,^o denote the non-relativistic field operators for the charged and for 
neutral pion fields, respectively. The propagator of the non-relativistic charged 
pion field is given by 



.(0|™,(.)*i(0)|0> = J ^ M. + pV2M.-p°-iO ■ 

The relativistic corrections due to the higher-order terms in Eq. (j4.2p are treated 
perturbatively. The non-relativistic propagator for the neutral pion is obtained 
by replacing M^r M^^o . 

The free non-relativistic field operators ^±,^q annihilate the vacuum. This 
property can be used to construct a theory that - from the beginning - conserves 
the number of pions. The interaction Lagrangian is then given by an infinite 
series of 4-pion local operators with an increasing number of space derivatives. 
In particular, in the 2-particle sector with zero total charge - spanned by the 
states |7r"''7r~) and Ivr^vr'^) - the interaction Lagrangian is written as 

£/ = ci^>;^$L^>+$- + C2(^>5.$L^>o^o + h-c.) + C3<I>J$J^>o$o + • • • , (4.4) 

where the ellipses stand for derivative terms. The vrvr scattering amplitude is 
calculated by using standard Feynman diagram techniques. To be specific, we 
consider the process 

7r+(pi)7r"(p2) vr+(p3)7r~(p4) . (4.5) 

Owing to the conservation of the number of pions, the structure of Feynman 
diagrams is particularly simple: to all orders, the pertinent Green function is 
determined by the bubble diagrams displayed in Fig. [3j In the CM frame P'^ = 
P1+P2 = (-P'^)O) , the contribution from Fig. (3]: is proportional to the loop 
integral 

dH 1 1 



J{F 



i(27r)4 + 12 /2M^ - PO + - iO + F /2M^ - P - iQ 

p, = y^MM-^M^)- (4.6) 
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Figure 3: Non-relativistic theory. Typical diagrams which contribute to the 
two-pion elastic scattering amplitude of the process ()4.5p . Solid (dashed) lines: 
charged (neutral) pions. Filled boxes and crosses denote derivative vertices and 
self-energy insertions, respectively. 

We have used dimensional regularization in intermediate steps of the calculation 
to tame ultraviolet divergences. Neutral pion loops are obtained with the replace- 
ment — > M^o. The contribution to the scattering amplitude for the process 
Eq. (j4.5p is obtained by putting Pq = 2 + = y^, where p denotes the 

pion three momentum in the CM frame. Therefore, the bubble graphs in Fig. [3] 
generate polynomials in the quantities 



The so constructed non-relativistic scattering amplitude reproduces the general 
low-energy structure of the relativistic amplitude in Eq. ()4.ip . The counterparts 
of the functions Ai{s), - ■ ■ , Di{s) are given in form of a power series in Pc,Pq, with 
coefficients that depend on the non-relativistic couplings ci , C2 , C3 , • • • . 

In order to ensure that the relativistic and the non-relativistic theories de- 
scribe the same physics at low energies, it remains to match the two theories, or, 
what is the same, to fix the non-relativistic coupling constants ci, C2, C3, • • • . The 
matching condition is formulated for the T-matrix elements. 



T^''"^{pi,P2;p3,P4) = [2Wa{pi) ■ ■ ■ 2Wd{pA)]^^^T'^R'^{pi,P2;P3,PA) , (4.8) 



where the subscripts R and NR label the relativistic and non-relativistic theories, 
and Wa{p) = \l M'^a + . The presence of the overall factor in the matching 
condition (j4.8p reflects the difference in the normalization of the one-particle 
states and the field operator in the non-relativistic and relativistic theories. It 
is understood that both sides of this equation are expanded in powers of the 
momenta pj. The matching should therefore be performed at a given order in 
the momentum expansion - it fixes the polynomial parts of the amplitudes in all 



Pc = [1 + 0{a^)\ , = Af^o [1 + 0{al)\ . 



(4.7) 
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physical channels. This is exactly the freedom one has in choosing the couplings 
of the non-relativistic Lagrangian. On the other hand, the non- analytic pieces 
proportional to cJc, ctq are reproduced automatically, according to analyticity and 
unitarity, which hold both in the relativistic and in the non-relativistic theories. 

This non-relativistic effective theory obeys power counting rules in a generic 
small 3-momentum p: Bubble diagrams with charged (neutral) pions running 
in the loop are proportional to pc (po)- Consequently, multi-loop diagrams are 
suppressed by pertinent powers of PcPo- It can be shown that the relativistic 
insertions and derivative couplings in the diagrams do not destroy the power 
counting. 

The diagrammatic expansion in this theory coincides with the effective range 
expansion. This can be seen most easily, if one assumes isospin symmetry = 
M^o. In this case, the bubbles vanish at threshold where Pc = Po = 0, and so 
do the contributions from the derivative vertices. This means that, to all orders, 
the threshold amplitudes are determined in terms of the non-derivative couplings 
ci , C2 , C3 . Using the matching condition ()4.8p , one can express these couplings 
through the vrvr scattering lengths with definite isospin, 

SM^ci = 47r(2ao + a2) + -- - , 
3M^C2 = 47r(a2 - ao) + ■ ■ ■ , 

3M^C3 = 27r(oo + 2a2) + -- - , (4.9) 

where the ellipses stand for isospin-breaking corrections. Analogously, the deriva- 
tive couplings in the Lagrangian can be expressed through effective ranges, shape 
parameters, etc. This property is ideally suited for describing hadronic atoms: 
the scattering lengths, which we want to extract from experimental data, turn 
out to be the parameters of the Lagrangian which will be used to describe the 
atoms. Consequently, the calculation of atomic observables in perturbation the- 
ory by using this Lagrangian will automatically generate a parametrization of the 
former directly in terms of scattering lengths. 

4.2 Including photons 

The inclusion of virtual photons in this framework is straightforward. First, one 
follows the paradigm of minimal coupling and replaces ordinary space-time deriva- 
tives of the charged pion fields by covariant ones. In addition, the Lagrangian 
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contains the kinetic term for free photons and a tower of gauge- and rotationaUy- 
invariant operators, which can be built from the electric E and magnetic B fields. 
For example, the kinetic term for charged pions becomes 

(j-\2 T-j4 DE ED \ 

+ + 8Mj =f ='"-6s^ +■ ^ 

where Dt^± = {dt^ieAo)^±,'D^± = V<I>-|- it ieA<I>-|- are covariant derivatives, e 
is the electric charge, and {Ao,A) denotes the photon field. Furthermore, hi is a 
new LEG, related to the electromagnetic radius of the pion, hi = M^{r'^) + 0(a). 

The power-counting at tree-level, which amounts to counting the number of 
3-momenta in a given Feynman diagram, can be carried out analogously to the 
case without photons. However, loop corrections in general lead to a breakdown 
of naive power-counting rules. This is a well-known problem, caused by the 
presence of a heavy scale in the Feynman integrals: loop integrals receive 
contributions from regions where the integration momenta are of the order of M^r, 
which cause a breakdown of the counting rules. On the other hand, the effect 
can be completely removed by simply changing the renormalization prescription 
in the non-relativistic EFT. This is so because the terms which break power 
counting behave like polynomials at low energy. Most straightforwardly, the goal 
can be achieved by modifying the prescription for the evaluation of Feynma n 
integrals. The pertinent modification is called "threshold expansion" (j99l . llOd ). 
A detailed descrip tion in the context of the hadronic atom problem can be found, 
e.g., in Refs. (|57l . Il09l ). In brief, the method boils down to Taylor-expanding 
the integrand in any Feynman integral in powers of the 3-momenta prior to 
performing the loop integrals in dimensional regularization. The expansion and 
the integration do not commute: it can be shown that the two results differ by just 
the above-mentioned polynomial contribution, which is absent in the threshold- 
expanded integral. Thus, applying threshold expansions to all loop integrals leads 
to a restoration of the naive power counting rules in the non-relativistic EFT. 

In principle, the matching condition in the presence of photons is again given 
by the relation Eq. (|4.8|) . On the other hand, the scattering amplitudes in the 
presence of real and virtual photons are infrared-divergent in perturbation theory. 
It is then natural to identify non-singular parts of the amplitude, which are more 
convenient for matching. At the accuracy needed here, it suffices to discuss 
the problem at order e^, for the charge-exchange process 7r"'"7r~ vr^vr*^. The 
structure of the scattering amplitude in the vicinity of the threshold |p| — > is 
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identical in the relativistic and in the non-relativistic case, 

e-^'-'^T-^-'^^ = ^ + e% In 1^ + T+-'00 + 0(p) , (4.11) 

I Pi -'"TT 

where p denotes the relative 3-momentum in the CM frame and Oc is the (infrared- 
divergent) Coulomb phase, 

The scale /x is generated by dimensional regularization, which is used to tame 
infrared and ultraviolet divergences. The coefficients 61^2 differ by a factor AM^ 
in the relativistic and in the non-relativistic theory. Finally, '^^ denotes 
the threshold amplitude, which is the counterpart of the scattering length in the 
presence of photons. It is infrared-finite. The matching condition at threshold 
reads 

7-+-;00^4^2^+-;00_ (4_^3^ 

Calculating the diagrams of the type shown in Fig. HI we arrive at an expression 
for Tj^^''^^ in terms of the non-relativistic couplings ci, C2, C3, • • • . At the order of 
accuracy we are working, only a finite number of diagrams contribute. The final 
result for the real part of the threshold amplitude is given by 

ReT+-°° = 2c2-C2ci^^ + ciC2^(^l-A(/i)-ln-fj +o(<^), 

(4.14) 

where A^r = — M^o, and A(/i) stands for the ultraviolet divergence originating 
from the diagram in Fig. UJi, 

A(/i) =/i2('i-3)|_l_ _p/(l) _1j^4^| _ (415) 

The ultraviolet divergence is removed in a standard manner, by renormalizing 
the coupling 02- 

The matching condition (j4.13p enables one to relate a particular combination 
of the couplings to the relativistic threshold amplitude At the accuracy 

needed here, higher-order terms in the momentum expansion of the amplitudes 
are not needed. 

Finally, we note that the non-relativistic couplings q contain both, strong and 
electromagnetic isospin-breaking corrections. According to the unified counting 
of the isospin-breaking effects, which was introduced in section [3l we write 

Ci = c, + Qcf ^ + (mrf - mu)^cf^ + o{6) , (4.16) 
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a b c d 




e f 

Figure 4: Representative set of diagrams contributing to the scattering amplitude 
of the process 7r"'"7r~ tt^t:^ at the order of accuracy we are working. The sohd 
and dashed hues denote charged and neutral pions, respectively, and the wiggled 
line denotes the Coulomb photon (the calculations are done in the Coulomb 
gauge). Transverse photons do not contribute at this order. Only the diagrams 
a,d,f contribute to the expression of the real part of the threshold amplitude, see 
Eq. dHH). 

where the bar denotes quantities taken in the isospin limit a = 0, rrid = rriu- The 
Ci can be related to scattering lengths and effective ranges in the isospin sym- 
metric world. On the other hand, the coefficients c-^'* , cp^ are fixed via matching 
to ChPT. 



4.3 Bound states 

The non-relativistic framework does not contain any new dynamical information 
about the behavior of the scattering amplitudes at low momenta, because it is 
constructed such that it reproduces the relativistic amplitudes. However, the 
non-relativistic framework is extremely useful when bound states are considered, 
because methods of standard quantum mechanics can be used to a large ex- 
tent. As all couplings in the non-relativistic Lagrangian have been fixed through 
matching of the scattering amplitudes, there are no additional free parameters 
left in the bound-state sector. Consequently, solving the bound-state problem 
in the non-relativistic theory, one can eventually express the observables of the 
bound states in terms of the parameters of the relativistic scattering amplitudes. 
We now describe the procedure. 

Hadronic atoms are shallow quasi-stable states formed predominately by the 
Coulomb force. In order to describe such states, it is convenient to use perturba- 
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tion theory, where the unperturbed solution corresponds to the purely Coulom- 
bic bound state. The full Hamiltonian of the system is constructed from the 
Lagrangian with standard methods. Here, we concentrate on that part of the 
Hamiltonian which is responsible for the next-to-Ieading-order term in the DGBT 
formula. Moreover, we confine f or si mplicity the calculation to the width of the 
ground state. As shown in Ref. (jl09l ). the pertinent Hamiltonian is 

H = Ho + He + H5 + Hk + • • • = Ho + He + V + • • • , (4.17) 

where tlo,c,S,R stand for the free non-relativistic Hamiltonian, Coulomb inter- 
action, strong interactions and the relativistic corrections to the pion kinetic 
energy. [It is convenient to use Coulomb gauge in the non-relativistic calcula- 
tions. In this gauge, transverse photons can be dropped completely, since they 
do not contribute at this accuracy to the width. The time-like photon field can 
be eliminated by using equations of motion, resulting into the static Coulomb 
potential acting between pions. Finally note that it is legitimate to use different 
gauges in the relativistic and non-relativistic theories, since only gauge- invariant 
quantities enter the matching condition.] The ellipses in Eq. (I4.17P stand for 
terms that do not contribute to the width of the ground state at next-to-leading 
order. The explicit expressions read 



H, 



= / d^x ^ cDt (x, 0) (u^a - $,(x, 0) 

a=±,0 ^ '^^ ^ 

He = /"d3x(i3y(^i(x,0)«I>_(x,0))^^(«&^(y,0)<&+(y,0)), 
4vr 7 |x - y| ^ 

= I d'x |-Ci ^>V$L$+«>„ - C2 (^^^L^*^ + h.C.) - C3 ($J)$o)'} , 

H^, = 1^3^ ^ ci>t(x,0)(^-^)$,(x,0). (4.18) 

It is seen that the Hamiltonian is amazingly simple: it contains three couplings 
ci,2,3 that need to be matched - all the rest is known. 

The pure Coulomb state is an eigenstate of the Hamiltonian Ho + H^. The 
resolvent Gc{z) = {z — Hq — He)~^ develops a tower of poles on the negative 
real axis in the complex z- plane, at z = En = 2M.,^ — a'^ M.,^ / An? , n = 1, 2, • • • (in 
the CM frame). The position of these poles coincide with the Coulomb binding 
energies. Once the perturbation V is switched on, the poles move from the real 
axis to the second Riemann sheet in the complex z-plane. The energy shift and 
width of a given state is defined by the real and imaginary parts of the shifted 
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pole position. Restricting ourselves to the ground state, we write 

r 

AE-i- = z-Ei, (4.19) 

The shift z — E-y of th e pole position can be consistently treated with the Feshbach 
formalism (jllTI . Ill8l ). A detailed dis c ussio n thereof in the context of hadronic 
atoms can be found in Refs. (j57l . 1 1041 . Il09l ). There, it is shown that the shift is 
given by the standard expression known from the Rayleigh-Schrodinger pertur- 
bation theory, 

z-E, = {^g\W + ^ Y1 '"^"i^l"' V + ---Wg), (4.20) 

where the sum over a runs over both the discrete and continuous spectra of the 
unperturbed Hamiltonian Hq + He, and denotes the ground-state vector. 

In momentum space, 

^gW= (k2+^2)2 ' 7 = l/rB = aMj2. (4.21) 

The center-of-mass (CM) motion is removed in the above matrix elements, which 
are then evaluated in the CM frame P = 0. 

It is instructive to first neglect the relativistic corrections H/j. At leading order 
in V, the shift of the pole position is real. Using Eqs. (|4.18p . ()4.20p and ()4.2ip . 
we get 

AS = -ii-l^ci + 0(V2), r = 0(v2). (4.22) 
Svr 

The matching condition displayed in Equation (j4.9p finally leads to 

AE = -la^M^{2ao + a2) + --- . (4.23) 



The ellipses denote contributions of order 6^. 

The decay width is of order . This leading term is generated by the contri- 
bution from the two neutral pion intermediate state in Eq. (|4.20p . The pertinent 
threshold is below the bound state energy - these states therefore generate an 
imaginary part in the energy shift. Since neutral pions do not feel the Coulomb 
potential, the sum over those intermediate states in Eq. (j4.20p merely yields the 
bubble integral with two neutral pions, similar to the one displayed in Eq. ()4.6p . 
The result for the width at this order reads 



r = -2Imz, 



(4.24) 
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where z is a solution to the equation 



4^ cUoiz). (4.25) 

Here, Jo{z) is given in Eq. (j4.6p . with M^r replaced by M^o. The equation ()4.25p 
has a solution on the second Riemann sheet only. The width becomes 

r = p ' C2 + -- - = - a p ' [ao-a2) +■■■ , (4.26) 

where p = 2M^o{Mt^ — M^o — M^a^/8). In the last step, the matching condition 
()4.9p was used. 

It is seen that the result for the width is of order 5^/^ at leading order. To 
work out the next-to- leading order terms, one has to include in Eq. (j4.20p con- 
tributions up to and including V'^, with V = H5 + Hr. According to power 
counting, the subsequent terms are suppressed by positive powers of 6. This can 
be seen e.g. from Eqs. (j4.6p and (j4.7p . showing that the charged and neutral bub- 
bles, evaluated at the bound-state energy = Ei, count as 0{5) and 0{5'^l^), 
respectively. This is a very important property of the non-relativistic EFT in 
dimensional regularization: at a given order in 6, only a finite number of terms 
in the perturbation series contribute. 

Finally, the result for the decay width of pionium up to and including terms of 
order (5^/^ reads 

where g{Ei) corresponds to the sum of diagrams where any number of Coulomb 
photons is exchanged between the charged pions. The explicit expression for this 
quantity is given by 

9{Ei) = ^ f 2 In a - 3 + K{p) + In ^ J • (4.28) 

Using the matching condition (j4.14p . one may finally express the decay width 
through the relativistic threshold amplitude of the process 7r"'"7r~ vr'^Tr'^, 

r = |gji(ReT+-°V(i + i^) + «(5^/^), 

K = ^(ao + 2a2)'-y (lna-l)(2ao + a2). (4.29) 

We note that the reference to the non-relativistic theory has completely disap- 
peared in the final result Eq. ()4.29p : the decay width is expressed through the 
relativistic threshold amplitude. 
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We have thus achieved the first main goal mentioned at the end of section 
[3j It remains to express the relativistic threshold amphtude in terms of isospin 
symmetric scattering lengths. Then, one can extract these from the measured 
lifetime of the pionium ground state. 



4.4 Scattering lengths 

The prediction for the scattering lengths 09.2 in Eq- (|2.3p concerns an isospin- 
symmetric paradise world - QCD at = rrid- In this world, there are no 
electromagnetic interactions. The light quark masses = m^, and the scale 
Aqcd are chosen such that = = 139.57 MeV, Mk = 493.68 MeV, 

F7r=92.4 MeV. The precise values of the heavy quark masses mc^b,t do not matter 
in the present context. On the other hand, the threshold amplitude, which occurs 
in the DGBT formula, concerns the real world, where rriu 7^ ma, a ^ 0. We are 
thus faced with the problem to relate that amplitude to the scattering lengths 
evaluated in the paradise world. 

The structure of the threshold amplitude at a 7^ 0, m„ 7^ md is 



— ReT+-^°° 
327r ^ 



ReT^ ' = ao — a2 + hi^rriii — rriu)'^ + + o{6) , (4.30) 



where the coefficients hi can be systematically calculated in the framework of 
ChPT. These calculations are carried out for the scattering amplitude, not for the 
bound state observables. Thus, the use of the non-relativistic approach enables 
one to separate bound state calculations from the chiral expansion. 

We outline the determination of hi^2 at order p^. The leading order Lagrangian 
of ChPT is 

^2 = — {d^Ud^U^ + 2BM{U + U^)) + CiQUQU^) , (4.31) 

where the unitary matrix U contains the pion fields, (• • • ) denotes the trace in 
flavor space, and 

M = diag{m^,md) , Q = | diag (2, -1) (4.32) 

are the quark mass matrix and the charge matrix, respectively. Finally, the 
constant C is related to the charged and neutral pion mass difference M^ — M^o = 
2e2C/F2. 

The threshold scattering amplitude n~^Tr~ — > tt^tt^ at order p'^ is given by 

+„;00_ s - Ml, _ 3M2 
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Prom this result, the expressions for hi at leading order can be read off, 

3A^ 



hi = 0{m) , /i2 



+ 0{m) . 
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(4.34) 



32a7rF2 

The d etails of the calculation at next-to-leading order can be found, e.g., in 
Refs. (jlOTI . Il09l ). The result for the width at next-to- leading order is 
2 



r = ^aV(ao - 02)^(1 + 5r) , Sr = (5.8 ± 1.2) x 10" 



(4.35) 



Note that the bulk of the total correction is generated by the leading-order 
term (|4.34p . which contains no free parameters. We expect that next-to-next-to- 
leading corrections will be completely negligible. Vacuun i pqlarizatiori has bee n 



investigated in pionium and/or other atoms in References (|57l . llOll . Il09l . Il 13 . 11201 ) . 

Using the scatt erin g lengths in Eq. ()2.3p . we arrive at the prediction for the 
pionium lifetime (|l09l ). 

1 



(2.9 ±0.1) X 10"^^ s. 



(4.36) 



The result of the ongoing measurement carried out by the DIRAC collaboration 
agrees with this value. 



r = {2mtoil) X 10-^^ s [DIRAC, Ref. jj] . 



(4.37) 



It is expected that the precision of the measurement improves in the near future, 
see Ref. (3). 



5 Pionic hydrogen and pionic deuterium 

The power and beauty of the non-relativistic effective Lagrangian approach is best 
demonstrated by the fact that the description of all hadronic atoms, which were 
mentioned in the introduction, proceeds very similarly to the pionium case just 
discussed. On the other hand, some of these bound systems are very different 
physically, and so are the results obtained. The important point is that the 
language used to describe these systems stays - with only minor modifications - 
always the same. 

As an example, we consider in this section the measurement of the 5-wave 
vrA^ scattering lengths aj^, in experiments on pionic hydrogen and pionic 
deuterium, which are performed by Pionic Hydrogen collaboration at PSI (|lol . 

. Measuring the energy shift and the width enables one to extract 
very accurate values of the real and imaginary parts of the elastic tt~p threshold 
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scattering amplitude, using pretty much the same technique as in the pionium 
case. Using unitarity and the measured Panofsky ratio finaUy allows one to to 
extract from data the real part of the threshold amplitudes for both, the elastic 
TT^p — > iT~p and the charge-exchange 7r~p — > Tr^n reactions. 

In the last step, the threshold amplitudes are again related to the pertinent 
scattering lengths in the isospin-symmetric world (cf. with subsection I4.4p . At 
leading order, the relation is 



TV p — >7r p 



T 



q-TT p- 



+ 



+ 



1 



4A, 



47r(l + Mjmp) \ 

1 f 9A^n 



ci 



(4/i + /2 



(5.1) 



167r(l + M^/mp) \mpF^ 

where nip denotes the nucleon mass, is the axial coupling constant of the 
nucleon and ci , /i , /2 are various (stro ng and e l ectro magnetic) LECs from the 



second-order pion- nucleon Lagrangian (jllOl . 11211 . |122| ) . (Following the tradition 
in the literature, we use the same notation for the LEG ci as in pionium. We 
hope that this is not confusing.) 

The difference between pionium and pionic hydrogen becomes visible by com- 
paring Eqs. (I4.30|4.34p and (jS.lh : whereas the pertinent isospin-correction for 
pionium at leading order is parameter- free, Eq. (15. ip contains the LECs ci, /i, /2, 
whose values are not established very well. The issue has been discussed in detail 
in Ref. (jSTI ) where, in particular, an update on the values of ci, /2 can be found. 
No reliable determination on the basis of experimental input is available for /i at 
present. This is the reason for a substantial uncertainty in the leading correction 
in the vrA^ case, which by far exceeds the experimental error in the measurement 
of the energy shift. 

Next, we turn to pionic deuterium, which allows one to extract the pion- 
deuteron threshold scattering amplitude. However, in this case the analysis is 
not yet complete: what one intends to finally obtain are the pion-nucleon scatter- 
ing lengths, which are related to the pion-deuteron amplitude through multiple- 
scattering theory. This is a very complicated issue, which has been extensively 



addressed in the past within the framework of potential mode 
lations in EE' 



M, m, m, 



67 



have been per formed as well (see, e.g., Refs. (|58l . l59l. iGOl. 



m, MM, 



s. Recently, ca. 



61 



71 



cu- 



72l)). This method allows one to largely reduce an 



uncontrolled systematic error in the resulting values of the vrA^ scattering lengths. 

The calculations within EET have shed new light on the importance of isospin- 
breaking corrections, a point which is obscure in potential models. Namely, the 
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pion-deuteron scattering length in the isospin hmit vanishes at leading order. For 
this reason, the isospin-breaking correction to this quantity, determined predomi- 
nately by short-range physics, turns out to be very large (6g). In the cont ext o f 
the pion-nucleon scattering, the same effect has been mentioned in Ref. (Il23l ). 
Further, the isospin-breaking correction to the pion-deuteron threshold ampli- 
tude contains the same virtually unknown LEG fi as the tt~p elastic scattering 
amplitude and is therefore determined with a large systematic error. 

The experiments on pionic hydrogen and pionic deuterium are complementary 
to each other. Namely, the data on pionic hydrogen alone determine the scattering 
lengths and aQ_,_ separately. The data on the energy shift of pionic deuterium 
provides an additional constraint on these two quantities. This can be seen by 
considering the bands in the {aQ^,aQ_^) plane, which correspond to the different 
observables. Measuring each of the following three observables: the energy shift 
and width of pionic hydrogen and the energy shift of pionic deuterium fix a 
particular combination of aQ_^ and Og^. Each combination corresponds to a band, 
whose width is determined by a combined experimental and theoretical error. If 
these three bands do not have a common intersection area in the {aQ_^, %+) plane, 
then either experiment or/and the theoretical interpretation of the data is not 
correct. 

The consistency check is made complicated by a large uncertainty pres ent in 
the values of LECs - most notably in ci and fi. Baru et al. in Ref. (|124| ) have 
proposed a procedure to partially circumvent the problem. The idea is based on 
the observation that the LECs ci and /i enter in the same combination in both, 
the pionic hydrogen and pionic deuterium energy shift (at leading chiral order). 
Introducing the quantity 

4tt{1 + M^/7np) V 

it is seen that ci and /i disappear at leading order from the expressions of the 
hadronic atom observables, if these are written in terms of a+ and aQ^_. Hence, 
various bands in the plot shown in Fig. [5] are much narrower than the pertinent 
bands that can be drawn in {a^j^, aQ_|_) plane (see Ref. (j57l )). In particular, Fig. [5] 
demonstrates that three bands still fail to pass this elaborate consistency test. 
However, it can be argued that next-order isospin-breaking corrections can be 
large and may change the above picture. In order to carry out a meaningful test, 
these corrections should be calculated at least to 0{p^) for all three observables. 
The possibility to perform the consistency check at higher accuracy, however, is 



«+ = 4+ + TZT^^rr^ZTT ( ^ ci - 2e2/i ) , (5.2) 
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Figure 5: The quantity a"*" defined in Eq. (|5.2p . plotted against a^^_. Three 
different bands emerge from the measurement of different observables. It is seen 
that with the isospin breaking corrections evaluated at 0{p^), the three bands 
1), 2) and 3) have no common intercept. 



not yet the end of the story. As it can be seen from Eq. (j5.2p , the relation of the 
isospin-symmetric scattering length a^^ to the quantity a_(_ does contain both, 
ci and fi. Thus, in order to determine aj^_ at a reasonable accuracy, one should 
find ways to estimate these LECs at the required precision. 

Finally, we mention that the extraction of the KN scattering lengths from 
the experimental data on kaonic hydrogen and kaonic deuterium 



H, Q, 0, 



18l . Iig . l2d ) bears many similarities to the pion-nucleon case. The analysis of 
the problem within three-flavor ChPT, however, is more complicated due to the 
large value of the strange quark mass. In addition, the presence of the sub- 
threshold A(1405) resonance leads to a large 5-wave scattering length. As a 
result, the deuteron problem can no more be treated purely perturbatively, and 
a partial re-summation of the multiple-scattering series should be considered. 
These very interesting issues, however, cannot be covered in the present review. 
The inte r ested rea der is referred to the original publications, e.g., Refs. (jSTI . IllSl . 

im, UM, I127I.I128I). 
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6 Summary points and future issues 
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i) Precise data on the energy levels and lifetimes of hadronic atoms enable one 
to extract various hadronic scattering lengths, provided that a systematic 
method to work out the relation between data and scattering lengths is 
available. 

ii) As we discussed in this review, a very convenient framework is provided 
by the non-relativistic effective Lagrangian approach. Its non-relativistic 
feature are used in intermediate steps only - at the end of the calculations, 
all observables are expressed in terms of the underlying relativistic theory, 
QCD+QED. 

iii) Despite the fact that various hadronic atoms observed in Nature are gov- 
erned by very different underlying physics, the same framework based on 
the non-relativistic effective Lagrangians applies - with minor modifications 
- to all of them. This is a beautiful demonstration of the potential and the 
flexibility of non-relativistic EFT. 

iv) To date, the conceptual problems of the general theory of hadronic atoms 
have been clarified to a large extent. Now, the focus shifts mainly to ap- 
plications. Among these, we mention the evaluation of a full set of isospin- 
breaking corrections at third order in pionic hydrogen and in pionic deu- 
terium. In addition, it would be a major breakthrough to present a sys- 
tematic calculation of the kaon-deuteron scattering length in terms of the 
threshold parameters of the KN interaction beyond the static approxima- 
tion. 
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